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ABSTRACT 



Context. Molecular data of extreme environments, such as Arp 220, but also NGC 253, show evidence for extremely high cosmic ray 

(CR) rates (10 3 - 10 4 x Milky Way) and mechanical heating from supernova driven turbulence. 

Aims. The consequences of high CR rates and mechanical heating on the chemistry in clouds are explored. 

Methods. PDR model predictions are made for low, n = 10 3 , and high, n = 10 5 5 cm -3 , density clouds using well-tested chemistry 
and radiation transfer codes. Column densities of relevant species are discussed, and special attention is given to water related species. 
Fluxes are shown for fine-structure lines of O, C + , C, and N + , and molecular lines of CO, HCN, HNC, and HCO + . A comparison is 
made to an X-ray dominated region model. 

Results. Fine- structure lines of [CII], [CI], and [OI] are remarkably similar for different mechanical heating and CR rates, when 
already exposed to large amounts of UV. HCN and H 2 abundances are boosted for very high mechanical heating rates, while ionized 
species are relatively unaffected. OH + and H 2 + are enhanced for very high CR rates f > 5 ■ 10~ 14 s -1 . A combination of OH + , OH, 
HiO + , HtO, and HaO + trace the CR rates, and are able to distinguish between enhanced cosmic rays and X-rays. 

Key words. Cosmic Rays - Mechanical Heating - Interstellar Medium - Excitation 
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Observations of molecular tracers have suggested that both 
star formation and an active galactic nucleus (AGN, an accret- 
ing supermassive black hole) can drive the physics and chem- 
istry of the central kpc of active galaxies (|S anders & Mirabel 
1 99 6: iGao et alJll997t iGao & Solomonll 19991: iBaan et aiT 2008 ; 
Ivan der Werf et al.ll2010h ~ Specifically, Ultra-Luminous Infrared 
Galaxies (ULIRGs), with infrared (8 - 1000 /mi) luminosities 
' £ir > 10 12 L , appear to contain large reservoirs of dense, 
n > 10 4 cm 4 , gas. This gas provides the fuel for (intense) star 
formation and black hole accretion. 

The irradiation by UV photons (6-13. 6 eV; creating so-called 
photo n-dominated regions, or PDRs, e.g.. lHollenbach & TielensI 
1999) from a starburst or by X-rays (1-100 keV, cre- 
ating so-called X-ray dominated regions, or XDRs, e.g., 
iLepp & Dalgar no 1996; [Malonev et al.ll 19961) from an AGN is 
expec ted to produce different sig n atures in the molecula r chem - 
istry dMeiierink & Spaansl 12005; Meij erink et al .1120061 120071) . 
Recently, mechanical heating due to supernova a nd stellar out- 
flow d riven turbulence has been explored as well (lLoenen et al.l 
2008), which shows that warm molecular gas may persist be- 
yond a few magnitudes of extinction in molecular clouds. 
Furthermore, the radiation feedback from ultraviolet photons, re- 
sulting in warm dust of more than 50 K, affect the thermodynam- 
ics of the molecular gas, possibly resulting in a top-heavy ini- 
tial m ass function (IMF; iKlessen et al.ll2007t iHocuk & Spaansl 
12011 . 

The analysis of the physical state of this irradiated gas allows 
one to constrain the evolutionary state of, ambient star formation 
rate in, and the importance of feedback for ULIRGs (Baa n et aD 



l20l0h . In this work, we focus on feedback that results from mas- 
sive stars. In particular, we consider UV emission from hot stars, 
the (elevated) production of cosmic rays from the resulting su- 
pernova remnants, and the heating that results from turbulent 
motions that are driven by travelling supernova blast waves. 

A very nice case in point for the feedback effects men- 
tioned above is Arp 220, the proto-typical ULIRG at 77 Mpc, 
which boasts a in frared luminosity of Lvr = 1.3 • 10 12 L 
(ISoifer et al.lfl9 87). Arp 220 is a merger and contains two nu- 
clei separated by about 370 pc in the east-west direction (e.g., 
iRodriguez-Rico et al.ll2005 |) . The nuclear star formation rate of 
Arp 220 is 50-100 M Q yr 1 (ISmith et alJll998h . 

Observations of Ar p 220 at a resolution of 0.3" (100 pc) by 
ISakamoto et all (120091) indicate P Cygni type profiles in HCO + 
4-3, 3-2 and CO 3-2. These authors interpret their results as a 
~ 100 km s _1 outflow that is driven by supernova explosions (or 
radiation pressure from the starburst). Interestingly, as many as 
49 supernova remnants (SNR s) have been detect ed with VLBI (1 
pc resolution) at 18 cm by (Lonsd ale et al.l20 06). confirming the 
formation and demise of many massive stars. The latter authors 
derive a supernova rate of 3 yr -1 for the western nucleus and 1 
yr 1 for the eastern nucleus, more than two orders of magnitude 
larger than the entire Milky Way. 

Limits on t he cosmic ray (CR) rate in our Milky Way are dis- 
cussed by, e.g JVan der Tak & van Dishoeckl(l2000l) . They found 
a value of £ = 2.6 + 0.8-10~ I7 s~' with upper limits that are 
3 to 5 times higher. In this paper we adopt ( = 5 ■ 10~ 17 s _1 
as the canonical Milky Way value. Highly elevated cosmic ray 
rates are seen in the rec ent Gamma ray observations of the star- 
burst galaxy NGC 253 (lAcero et all2009h with the High Energy 
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Stereoscopic System (H.E.S.S.) array. The observed gamma ray 
flux indicates a cosmic ray density that is three orders of magni- 
tude larg er than our Milk y Way. 

Papadopoulos (2010) suggests that ULIRGs generally have 
cosmic ray densities Uqr ~ 10 3 - 10 4 t/ cli MW, and discusses 
the impact on the ionization fraction and thermal state of the 
gas, and the consequences for the initial conditions for star- 
formation. In this paper, we first discuss the detailed effects 
of enhanced cosmic ray rates on both the thermal balance and 
chemistry (Section 0. Then, we will highlight the global con- 
sequences on the chemistry and suggest potential diagnostic 
species (Section[3]l, which is followed by the observational con- 
sequences (Section HJ. In Section [3] we will conclude with the 
key diagnostics needed to derive quantitive and qualitative prop- 
erties of interstellar medium (ISM) exposed to large amounts of 
cosmic rays, and show how it differs from gas exposed X-rays. 



2. PDR models 

We have constructed a set of PD R models (se e Tab le [TJ 
from the code as desc ribed by iMeiierink & Spaani d2005l) and 
Meij erink et al.l d2007l) . We consider a low (n — 10 3 cm -3 ), and 
a high density (n = 10 5 ' 5 cm" 3 ) phase. The low density phase 
of the ISM is generally not directly connected to starformation. 
Therefore, the adopted incident UV field in the low density mod- 
els (Go = 10 3 ) is chosen a factor 100 lower than the high density 
phase (Go = 10 5 ) in order to reflect this and essen t ially c orre- 
sponding to Model 1 and 4 in Meijeri nk & Spaansl (120051) . We 
varied the incident cosmic ray rates between £ = 5 • 10~ 17 s" 1 
and £ = 5 • 10~ 13 s _1 . The highest cosmic ray rate considered 
here is two orders of magnitude higher than in Meij erink et al.l 
(2006). We adopt a Solar metallicity Z = Z Q and an abundance 
ra tio [Cl/TOl = 0.4, wh ich is based on the average abundances 
of lAsplund etall (|2005) and Jenkins (2004), and summarized in 
Table 2 of Meijerink & Spaans l2005l 

We also consider the effect of mechanical heating in the 
high density case, specifically r mec h = 3.0 • 10~ 18 and 2.0 • 
10~ 17 erg crrT 3 s _1 , which correspond to star formation rates of 
about S F R = 140 and 950 M yr , respectively, for a Salpeter 
IMF (see lLoenen et alj |2008, for details). Mechanical heating is 
only considered in the high density case, n > 10 4 - 10 5 cm -3 , 
since we expect this dense phase to be directly responsible for 
star formation and thus be affected directly by the supernovae 
from massive stars. We assume that the low density phase is 
only affected by cosmic rays. These are able to travel through 
large columns (/Vh ~ 10 24 cm" 2 ), since they have small absorp- 
tion cross sections, and are likely to impact the ISM on larger 
scales. 

The effects on the thermal balance and chemistry are dis- 
cussed for Model sets la (Go = 10 5 ; n = 10 55 ctrT 3 ), lb 
and lc (same as Model set la, with F mec h = 3.0 • 10~ 18 and 
2.0 • 10" 17 erg ctrT 3 s , respectively) and 2 (Go = 10 3 ; n = 
10 3 cm' 3 ) in order to highlight the differences between the high 
and low density phases of the ISM. 

2.1. Thermal balance 

The surface temperature for the high density model without 
boosting the CR rate (i.e. f = 5 • 10~ 17 s) is higher than for 
the low density model. This may seem a little counterintuitive, 
since Go/«=l (0.3) for the low (high) density model, and as 
a result there is more energy available per particle in the low 
density model. Besides that, cooling is more efficient at high 



Table 1. Model parameters 



Metallicity [Z Q ] 


1.0 








ratio 


n 1 








CR rates [s 1 ] 


5- 10- 17 , 


5 ■ 10~ 16 , 5 • 10-' 5 , 


5 • 10" 14 , 5 • 


io- 13 


Model set 


Density 


Radiation field 


Tmech 






[cirT 3 ] 


[Go] 


[erg cirT 3 




la 


10 5i 


10* 


0.0 




lb 


10 55 


10 3 


3.0- 10' 


-18 


lc 


10 55 


10 s 


2.0 • 10' 


-17 


2 


10 3 


10 3 


0.0 





densities. However, at high radiation fields (Go > 10 , and 
n > 10 4 - 10 5 cm" 3 ), the gas cooling thermalizes (A oc n in- 
stead of n 2 ), while photo-electric heating becomes more effec- 
tive because of the higher electron density, and more effe ctive 
than F ~ n ("e.g. JKaufman et alJfl999t IMeiierink et al.ll2007l) . 

The sensitivity of the thermal balance to cosmic rays is dif- 
ferent for the high and low density phase (see Fig.Q]). In Model 
set la, only very high CR rates (£ > 5 ■ 10~ 14 s~') are able to 
increase the temperature significantly (~ 25 to 200 percent) in 
both the radical region (low Ay) as well as the shielded molec- 
ular region. In Model sets lb and lc, mechanical heating in the 
shielded region of the cloud dominates over cosmic ray heating. 
Here, the temperatures overlap for all cosmic ray rates, reach- 
ing values T ~ 100 - 200 K for Model set lb and T ~ 400 K 
for Model set lc. So even for moderate mechanical heating rates 
(in this case corresponding to a SFR=140 M /yr), the CR rate 
has no additional effect on the thermal balance in the shielded 
regions of clouds. In Model set 2, temperatures are increased for 
all CR rates at all column densities, up to an order of magnitude 
with the highest adopted CR rate. The effect on the thermal bal- 
ance is larger, because cooling is simply less efficient at lower 
densities (A oc « 2 ), compared to the heating rate that is propor- 
tional to the density, (T oc n). This makes the low density mod- 
els more sensitive to a larger heating rate. This is enhanced by 
the fact that the fractional electron abundance is increased more 
easily as well at lower densities (k rec oc n 2 ), r esulting in a higher 
heating efficiency (cf., Bake s & Tielensl[l994l) . and higher heat- 
ing rates. 

2.2. Chemical balance 

Electron abundances (Fig. |2): The effect of cosmic rays is less 
apparent in the unshielded region (tVh 55 3 ■ 10 21 ctrT 2 ) than in the 
molecular region (/Vh Z 10 22 ctrT 2 ) of the cloud, since charged 
species in the unattenuated part of cloud are also produced 
through photoionization by the stellar radiation field. Here, the 
fractional abundance (which is defined as x\ - «i/n(H + 2H2) of 
electrons varies between x e ~ 10~ 4 - 10~ 3 (10~ 4 - 10~ 2 ) in the 
high (low) density models, while the CR rates vary over 4 or- 
ders of magnitude. In the molecular region, however, an order of 
magnitude higher CR rate gives a similar rise in the electron frac- 
tional abundance, giving ranges x e ~ 10~ 9 - 10~ 5 and 10~ 7 - 10" 2 
for the high and low density models, respectively. Adding me- 
chanical heating (lb and lc) does not change the electron abun- 
dance much. The difference in fractional electron abundance be- 
tween the highest and lowest CR rate is somewhat larger than for 
the models without mechanical heating. 

H and Hi abundances (Fig. |2): Atomic hydrogen is by 
far the dominant species in the radical region (/Vh < 2 - 3 • 
10 22 cm -2 ), but when the attenuating column becomes larger, 
molecular hydrogen becomes more and more abundant. For the 
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Fig. 1. Temperatures for all models sets (see Table [TJ for cosmic ray rates £ = 5 • 10 17 (solid line), 5 • 10 16 (dotted), 5-10 15 
(dashed), 5 • 10~ 14 (dot-dashed), and 5 • 10~ 13 s _1 (dotted-dashed). 
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Fig. 2. Abundances of H, H2, e , C + , C and CO for all models sets (see TableQ]). Line styles corresponding to the different CR rates 
are the same as in Fig.Q] 



canonical CR rate (( = 5 ■ 10 s ), the atomic hydrogen abun- 
dance drops below xh ~ 10~ 4 (10~ 2 ) for the high (low) den- 
sity models. The H to H2 transition becomes less and less com- 
plete for higher CR rates, which is due to CR reactions such as 
H 2 + CR -» H + H + + e and H 2 + CR -> H+ + e . Atomic 



and molecular hydrogen are of similar abundance for the highest 
adopted CR rate in the high density case, while atomic hydrogen 
is dominant for CR rates £ > 5 • 10~ 15 s _1 in the low density case. 

C + , C and CO abundances (Fig. |2): Similar trends are 
seen for the C + /C/CO transition. For the lowest CR rate = 



4 



R. Meijerink et al.: Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating. 



5 • lfr 17 s- 1 ), there are three well defined zones where one 
species is most abundant. Increasing CR rates result in a larger 
neutral carbon abundance in the unattenuated part of the cloud 
by a factor of ~ 3 (~ 30) in the high (low) density case. However, 
the abundances are still low compared to C + . Again most ef- 
fect is seen at large column density (A^h > 10 22 cm -2 ). The C + 
and C abundances increase very fast with the CR rate. In the 
high density case, the C + and C abundances become 10~ 6 - 10~ 5 
and 3 - 6 • 10~ 5 , respectively, depending on the amount of addi- 
tional mechanical heating. There is no full transition to CO for 
£ > 10~ 15 s in the low density case. Neutral carbon is the dom- 
inant species at high column density f = 5 • 10~ 15 s _1 , and C + 
has become the dominant species at f > 5 • 10~ 14 s , and at the 
same time there is almost no CO left (xco < 10~ 6 ). 

H + , Hi, and //J abundances (Fig. The main drivers of 
chemistry in PDRs are ionic species, starting with species H + , 
HJ, and HJ. H + and HJ are mainly (directly or indirectly) pro- 
duced by cosmic ray ionization (H + CR — » H + + e ; H2 + CR — > 
HJ + e~ or H + + H + e or H + H). UV photons of energy larger 
than 13.6 eV are all readily absorbed in a very small column, 
forming an HII region in front of the PDR. H2 can only be pho- 
todissociated to two hydrogen atoms (with an UV excited W 2 as 
intermediate), no ionic species are produced by direct photoion- 
ization of H2. Therefore, large effects on abundances of ionic 
species are seen by varying the cosmic ray rate. 

In the unshielded region, the abundance increase of H + and 
HJ scale with the ionization rate. H + is directly produced by 
cosmic ray ionization. The HJ abundance is balanced by the 
charge exchange reactions H + + H2 <-> HJ + H. The reaction 
to the right has an energy barrier of T ~ 20000 K, and is only 
efficient at temperatures of a few thousand Kelvin. Therefore, 
the HJ abundance is still able to increase with CR rate: the 
smaller H2 abundance at higher CR rates is counteracted by a 
temperature increase (see Fig.Q]). At slightly larger column den- 
sities (A^h ~ 3 • 10 21 cm -2 ), the temperature drops and only 
HJ + H — > H + + H2 is efficient, while simultaneously the H 
to H2 transition occurs. This together results in an abundance 
maximum of H + , and also explains why the HJ abundance is not 
increasing as fast as the H2 abundance. 

After the H to H 2 transition (Nh > 5 • 10 21 cm -2 ), the over- 
all ionization degree drops, which coincides with a decrease in 
temperature. This is mostly explained by recombination reac- 
tions (e.g., H + + e — > H + Photon), that are more efficient at 
lower temperatures. In this region, the increase of H + with cos- 
mic ray rate is faster than linear. The H to H2 transition is less 
complete for higher CR rates (Fig. 13, giving a larger reservoir 
of atomic hydrogen, that can potentially be ionized. 

HJ is formed and destroyed by HJ + H2 — * HJ + e~ and 
Hj + e~H — > 3H and H2 + H, respectively. At low column den- 
sities, the HJ abundance is dropping with CR rates due to en- 
hanced electron and reduced H2 abundances. At higher column 
densities, HJ increases for large CR rates, but not one to one 
with the CR rate, due to more efficient recombination and lower 
H2 abundances (especially in the low density models). 

CH + , C7/J, and CH abundances (Fig. In the unshielded 
region of the cloud, the formation of CH + occurs mainly through 
C + + Hj — * CH + + H, where H* 2 is vibrationally excited molec- 
ular hydrogen. Therefore, the CH + abundances do not change 
much with cosmic ray rate in the unshielded region of the cloud. 
However, at larger column densities (A^h ~ 10 22 cm' 2 ), the 
main route is through HJ + C. At high density, both the HJ 
and C abundance increase orders of magnitude for the larger 
CR rates, and as a result the CH + abundance as well, despite 
the increased electron abundance. In the low density case, the 



CH + abundances are very similar for all CR rates because there 
is no boost now through C + HJ. CHJ formation occurs mainly 
through CH + + H2 — > CHJ + H and follows the abundance of 
CH + closely, except where molecular hydrogen is not the domi- 
nant species, i.e., the radical region in the high density case and 
for very high CR rates in the low density case. 

CH can be formed by the neutral-neutral reactions C + H2 — > 
CH + H and C + OH -> CH + O, which have energy bar- 
riers of T ~ 11700 and 14800 K, respectively, and the re- 
combination reaction CHJ + e~ — > CH + H. Destruction hap- 
pens through UV photodissociation and photoionization at small 
column densities, and by fast charge exchange reactions (e.g., 
H + + CH -> CH + + H) at all column densities. At small col- 
umn densities (Nh 5 10 21 cm -2 ), the CH abundances are fairly 
similar, and UV photoionization dominates the destruction, and 
neutral-neutral reactions formation. The small differences can be 
explained by small deviations in H2 and C abundances. At col- 
umn densities A^h ^ 5 • 10 21 cm -2 , the formation is dominated by 
recombination of CHJ, and the abundance is mainly determined 
by the CR rate. 

Water chemistry (Fig. H2O and related species, OH + , 
OH, H20 + , and H30 + are significantly affected by cosmic rays. 
First we consider H2O. It can be formed through the neutral- 
neutral reaction H2 + OH — > H2O + H, for which temperatures 
T > 200 - 300 K are required. Cosmic rays are not able to in- 
crease the temperature to these temperatures, not even the high- 
est cosmic ray rates. Another important route, especially when 
gas is significantly ionized by enhanced cosmic ray rates, is, e.g., 
HJ + O -> OH + + H, then OH + + H 2 -> H 2 + + H, H 2 + + H 2 -> 
H 3 + + H, followed by recombination, H 3 + + e~ -> H 2 + H. 
Destruction is possible through UV photodissociation, charge 
exchange, proton transfer, ion-molecule reactions and cosmic 
rays. Important routes are He + + H2O — > OH + He + H + and 
H + + H2O — > H20 + + H, and varying the cosmic ray rates af- 
fect these abundances. Therefore, at intermediate column den- 
sities Nh ~ 3 • 10 21 to 10 22 cm~ 2 in model sets la and lb, the 
water abundances increase for higher CR rates. In these mod- 
els, the temperatures are not high enough to form H2O through 
neutral-neutral reactions. When the temperature is high enough 
to form water (model set lc) through the neutral-neutral reac- 
tions, a higher ionization degree is a drawback (H2O is destroyed 
by, e.g., He + ), and the water abundance decreases for very high 
CR rates > 5 • 10~ 14 s ). In model set 2, there is not enough 
H2 at the higher cosmic ray rates to form water through ion- 
molecule reactions. 

OH is formed through H2 + O — > OH + H or destruction 
of water (see previous paragraph). In the radical region where 
temperatures are signficant, formation through neutral reactions 
dominates. In the high density models they are closely connected 
to the H2 abundance, and here the OH abundance is slightly 
lower for the high CR rates. In the low density case, the OH 
abundance increases with CR rate in the radical region, since the 
temperature increases about an order of magnitude. When tem- 
peratures are lower, i.e., Nh <: 5 • 10 21 crrT 2 , the contribution of 
water destruction becomes more important. OH abundances are 
significantly larger in model set la and lb, when CR rates are 
higher. In model set lc, formation through neutral-neutral reac- 
tions becomes more important, and the effect of an enhanced CR 
rate less evident. 

OH + is formed by H + + OH -> OH + + H and H + + O 
+ + H followed by 6 + + H 2 -> OH + + H or by photoionization 
of OH. At small column densities, where OH + abundances are 
very similar for different CR rates, the abundance is determined 
by the UV ionization of OH. At large column densities (A^h S 
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5 • 10 21 cirT 2 ), the UV is mostly gone, and here the H + driven 
formation is more important, which is enhanced by cosmic rays. 
H20 + closely follows the OH + abundance at these high column 
densities. It is mainly formed through H2 + OH + — » H20 + + H, 
and also very much enhanced at higher CR rates. 

Nitrogen chemistry (Fig. |5); N + has an ionization potential 
of 14.53 eV, which is too high for the available UV photons. 
This species shows orders of magnitude variations for all model 
sets. The species is strongly affected by the increased cosmic 
ray rate. NH + and NHj are closely connected to the cosmic ray 
rate at high column densities, because the formation is driven by 
N + H+ and N + + H 2 . 

HCN and HNC do not show a strong response to enhanced 
CR rates in the high density case without mechanical heating, 
and the abundances remain similar for all CR rates. Adding me- 
chanical heating (lb), however, drives the HNC + H -> HCN + H 
reaction that has a temperature barrier of T ~ 200 K, and HCN 
has an abundance that is two orders of magnitude higher than 
HNC. In the low density case (2) both the HCN and HNC abun- 
dances, which are similar, decrease fast with increasing cosmic 
ray rates. This is a direct results of the decrease of H2 and H^ , 
which are important in the formation of these molecules. 

HCO + andCS (Fig.®: For Model sets la and lb, the HCO + 
abundances are enhanced with increasing CR rates. The forma- 
tion occurs through H+ + CO -> HCO + + H 2 , which is fast 
at all temperatures, and is aided by CH + O — > HCO + + e~ at 



high temperatures. Destruction happens through recombination 
with electrons but also through reactions with molecules such as 
H 2 0, OH, and HCN. In Model set lb, the abundances of these 
species are especially enhanced for the lower ( f = 5 • 10~ 17 
and 5 • 10~ 16 s~') cosmic ray rates as discussed in the previ- 
ous paragraphs, and suppress the HCO + abundance compared to 
the models without mechanical heating. In the low density mod- 
els (2), where the full transition to H2 does not occur, the H3 
abundance is automatically suppressed as well, giving smaller 
abundances with increasing cosmic ray rates. 

The sulfur chemistry contains many reactions with activa- 
tion barriers. CS has two formation channels, the first through 
S + CH -> CS + H and the second through S + H 2 -» HS +H (with 
an energy barrier of T ~ 6000 K) followed by HS+C -> CS + H. 
CS is most abundant at column densities Nh <; 5 • 10 21 cm" 2 , 
where temperatures are not such that the second channel is very 
effective, although not insignificant in Model set 2. The forma- 
tion of CH on the other hand is dependent on Ht , of which the 
abundance is highly dependent on CR rate. Destruction hap- 
pens through ion-molecule reactions with, e.g, He + , H + , H^, 
that are also dependent on CR rate. The abundance is a very 
complex interplay between temperature a nd ionization rates (cf., 
iLeen & Graft! 1 9881: iMeiierink et alJ2008h . The abundance of CS 
increases with CR rate in Model set la and lb, while Model set 
lc shows a drop. Molecule formation becomes more and more 
ineffective in the low density case, and the CS abundance de- 
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Fig. 4. Abundances of OH + , H20 + , H30 + , + , OH, and H2O for all models sets (see Table [T). Line styles corresponding to the 
different CR rates are the same as in Fig.Q] 



creases significantly with CR rates, similar to, e.g., HCN and 
HNC. 



3. Integrated column densities 

3.1. Trends 

The integrated column densities for the total cloud (A^h = 3 • 
10 22 cm~ 2 ) are shown in Tables[2]through[5]for the species dis- 
cussed in the previous section. The main results are summarized 
in Fig. [7] 

Hydrogen chemistry (top row Fig. \7}: The integrated col- 
umn densities of atomic and molecular hydrogen vary over a fac- 
tor of 3 (10) for the high (low) density models over the whole 
range of cosmic ray rates. The ionic species of hydrogen, H + , 
Hj, and H3 show much larger variations. The column densities 
of N(H + ) vary from ~ 10 16 -10 18 ctrT 2 (10 17 - 10 20 ctrT 2 ) for the 
high (low density models). H^ increases (decreases) 4 (2) orders 
of magnitude in the high (low) density case. 

Water chemistry (row 2 and 3 of Fig. \7}: The + column 
densities show a nice correlation with CR rate, and increases 2 
to 3 orders of magnitude over the range of considered CR rates. 
Although there are large variations in absolute column densities 
for OH and H2O between the different models, the OH/H2O ratio 
increases with cosmic ray rate, from N(OH)/N(H20) ~ 4 • 10~ 3 
(~ 0.2) to ~ 1 - 3 (~ 400) for the high (low) density models. 



The OH + and H20 + show a strong response to the CR rate 
for £ > 5 ■ 10~ 15 s _1 in the high density case. Above this thresh- 
old, the column density increases with a factor ~ 10-20 for 
OH + and ~ 50 - 60 for H20 + . In the low density case, the col- 
umn densities for both species show a steep rise with CR rate 
up to ( ~ 5 ■ 10~ 14 s _1 , after which it levels off for OH + and 
decreases for H20 + . The trends of HsO + generally show an in- 
crease (decrease) with CR rate in the high (low) density case. 
The response is more direct when high mechanical heating is 
added and temperatures are higher. 

CO + and HCO + (row 4 in Fig. 0: HCO + does not show 
an obvious trend with increasing CR rates in the high density 
case. The HCO + abundance anti-correlates with CR rate in the 
low density case. It does not seem to be a tracer that can be 
easily used as a CR tracer, but it does provide information when 
used in combination with e.g. the water chemistry. CO + is very 
irresponsive to the CR rates at high density except for very high 
CR rates (f > 5 ■ 10~ 14 s _1 ). In the low density case, N(CO + ) 
does not show a very obvious trend. 

C + , N + , and C (bottom row in Fig. [7}: Neutral and singly 
ionized carbon have very constant column densities over the 
range of CR rates considered. N + on the other hand correlates 
very well with the CR rates, and therefore the N + /C or N + /C + 
column density ratio can be diagnostic in estimating the CR rate. 

Other species: Despite the fact that many species show large 
local variations in a cloud when varying CR rates, it does not al- 
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Column density [cm 2 ] 



Fig. 5. Abundances of N + , NH + , NH^, HCN, and HNC for all models sets (see Table[TJ. Line styles corresponding to the different 
CR rates are the same as in Fig.Q] 




Column density [cm 2 ] 



Fig. 6. Abundances of HCO + and CS for all model sets (see Table [TJ. Line styles corresponding to the different CR rates are the 
same as in Fig.Q] 



ways result in significant changes in the integrated column den- 
sities in high density clouds. Examples are CH + and CHj which 
are the same within 5 percent over the whole range. In the low 
density models, variations are larger, but the integrated column 
densities are 1 or 2 orders of magnitude smaller. There are also 



species that are not easy to interpret, such as NFL and NFFJ\ The 
trends show a minimum or a maximum, when varying the CR 
rates. A number of species we considered have a strong depen- 
dence on temperature and ionization rate. Especially CS exhibits 
are very complex interplay, and does not show a particular trend. 
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Fig. 7. Column densities of diagnostic species. 



3.2. Uncertainties 

Molecular hydrogen, H2, is a key species in the chemistry of 
the ISM. The formation of H2, which mainly occurs on grain 
surfaces at these densities and metallicity, depend on the spe- 
cific properties and temperatures of the grains, and is s ubject 
to a considerable uncertainty (cf., Cazaux & Tielens 2004). As a 
result the H2 abundance is uncertain as well, especially when the 
gas is not predominantly molecular, i.e., in the radical region. 

Therefore, the column densities of key species such as + , 
OH, OH + , and H20 + might be affected, especially when the in- 



tegrated column density is dominated by the abundances in the 
radical region. In order to investigate this, we lowered the H2 
formation rate by a factor 3, and here we summarize the results 
for some of the species, considering Model set la: 

H2, //J, and H + : The abundance of H2 is lower by a factor of 
~ 3 in the unshielded region of the cloud, and the H to H2 tran- 
sition occurs now more gradual. In the molecular region, there 
is only a noticeable effect when there is no full transition to H2, 
which is for the highest considered CR rate, f = 5 • 10~ 13 s _1 . 
The column density of atomic hydrogen is ~ 10 percent higher 
for the lowest up to ~ 80 percent higher for the highest CR rate. 
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Species Cosmic ray rate (s ') 





5(-17) 


5(-16) 


5(-15) 


5(-14) 


5(-13) 




Column densities (cm" 






H 


4.2(21) 


4.2(21) 


4.3(21) 


5.2(21) 


1.2(22) 


H 2 


1.4(22) 


1.4(22) 


1.4(22) 


1.4(22) 


1.0(22) 


H + 


4.9(15) 


7.0(15) 


2.7(16) 


1.9(17) 


1.3(18) 


H 2 + 


2.7(09) 


2.1(10) 


2.0(11) 


2.0(12) 


1.6(13) 


H+ 


4.1(12) 


4.3(13) 


4.1(14) 


3.2(15) 


3.9(15) 


e 


8.1(17) 


8.1(17) 


8.4(17) 


1.0(18) 


2.5(18) 


CH + 


2.8(12) 


2.7(12) 


2.7(12) 


2.4(12) 


2.0(12) 


CH+ 


8.7(11) 


8.7(11) 


8.7(11) 


8.4(11) 


9.4(11) 


CH 


3.1(13) 


3.3(13) 


5.2(13) 


1.1(14) 


8.7(13) 





6.6(18) 


6.4(18) 


6.4(18) 


6.8(18) 


9.3(18) 


+ 


1.4(12) 


2.1(12) 


8.2(12) 


5.7(13) 


3.5(14) 


OH + 


1.6(12) 


1.6(12) 


2.1(12) 


6.3(12) 


3.4(13) 


H 2 + 


5.1(11) 


5.4(11) 


8.0(11) 


3.3(12) 


2.8(13) 


H 3 0+ 


1.6(14) 


3.1(14) 


4.2(14) 


6.1(14) 


1.7(14) 


OH 


5.7(14) 


1.5(15) 


8.1(15) 


4.8(16) 


6.7(16) 


H 2 


1.4(17) 


2.9(16) 


5.4(16) 


9.1(16) 


2.0(16) 


C + 


7.6(17) 


7.6(17) 


7.6(17) 


7.7(17) 


1.0(18) 


C 


4.7(17) 


5.0(17) 


5.4(17) 


6.3(17) 


1.8(18) 


CO 


3.2(18) 


3.2(18) 


3.2(18) 


3.1(18) 


1.7(18) 


CO + 


5.7(10) 


5.7(10) 


5.9(10) 


8.3(10) 


7.9(11) 


HCO + 


5.9(13) 


6.7(14) 


1.5(15) 


2.3(15) 


2.3(14) 


N + 


4.0(10) 


4.0(11) 


4.3(12) 


5.0(13) 


3.9(14) 


NH + 


3.5(11) 


3.5(11) 


3.3(11) 


2.9(11) 


3.9(11) 


NH + 


1.3(11) 


1.3(11) 


1.4(11) 


2.3(11) 


5.7(11) 


HCN 


4.3(13) 


5.7(13) 


5.6(13) 


1.2(14) 


3.2(13) 


HNC 


2.9(13) 


3.5(13) 


7.6(13) 


1.8(14) 


5.0(13) 


CS 


3.1(13) 


8.4(13) 


2.5(14) 


1.6(14) 


3.9(12) 



Note - Numbers in parentheses: a(b) = a ■ 10". 



Table 3. Column densities for Model set lb 



Species Cosmic ray rate (s ') 





5(-17) 


5(-16) 


5(-15) 


5(-14) 


5(-13) 




Column densities (cm" 






H 


5.1(21) 


5.1(21) 


5.3(21) 


6.2(21) 


1.4(22) 


H 2 


1.4(22) 


1.4(22) 


1.4(22) 


1.3(22) 


9.3(21) 


H + 


1.2(16) 


1.5(16) 


3.8(16) 


2.2(17) 


1.4(18) 


H+ 


2.7(09) 


2.0(10) 


2.0(11) 


1.9(12) 


1.5(13) 


H+ 


1.5(12) 


1.5(13) 


1.5(14) 


1.9(15) 


8.6(15) 


e 


7.8(17) 


7.8(17) 


8.1(17) 


1.0(18) 


2.3(18) 


CH+ 


5.7(12) 


5.7(12) 


5.6(12) 


5.5(12) 


5.2(12) 


CH+ 


1.9(12) 


1.9(12) 


1.9(12) 


1.9(12) 


2.1(12) 


CH 


2.9(13) 


3.8(13) 


4.5(13) 


7.4(13) 


6.2(13) 


O 


3.2(18) 


3.2(18) 


3.4(18) 


5.1(18) 


7.8(18) 


o + 


3.7(12) 


4.4(12) 


1.1(13) 


6.6(13) 


3.9(14) 


OH + 


5.2(12) 


5.3(12) 


5.8(12) 


1.1(13) 


5.9(13) 


H 2 + 


1.6(12) 


1.6(12) 


1.8(12) 


4.2(12) 


6.1(13) 


H,0 + 


3.3(13) 


1.5(14) 


8.5(14) 


3.2(15) 


1.8(15) 


OH 


3.1(16) 


3.2(16) 


3.8(16) 


8.5(16) 


2.1(17) 


H 2 


5.3(18) 


5.1(18) 


4.5(18) 


2.3(18) 


3.4(17) 


C + 


7.2(17) 


7.2(17) 


7.2(17) 


7.4(17) 


8.5(17) 


C 


1.2(17) 


1.8(17) 


3.4(17) 


5.3(17) 


1.1(18) 


CO 


2.3(18) 


2.5(18) 


2.9(18) 


3.1(18) 


2.5(18) 


CO + 


2.4(11) 


2.4(11) 


2.4(11) 


2.5(11) 


9.2(11) 


HCO+ 


1.9(13) 


3.9(13) 


2.2(14) 


1.9(15) 


6.4(15) 


N + 


6.8(10) 


6.8(11) 


6.9(12) 


6.6(13) 


4.4(14) 


NH + 


1.4(11) 


1.3(11) 


1.3(11) 


1.2(11) 


2.7(11) 


NHJ 


3.5(10) 


3.5(10) 


4.1(10) 


1.0(11) 


6.1(11) 


HCN 


8.0(17) 


8.2(17) 


4.4(17) 


9.5(16) 


3.5(15) 


HNC 


1.1(16) 


8.0(15) 


5.5(15) 


1.2(15) 


8.0(13) 


CS 


4.6(14) 


3.2(14) 


6.3(13) 


3.5(13) 


8.9(12) 



Species Cosmic ray rate (s ') 





5(-17) 


5(-16) 


5(-15) 


5(-14) 


5(-13) 




Column densities (cm" 






H 


4.2(21) 


4.3(21) 


4.3(21) 


5.1(21) 


1.3(22) 


H 2 


1.4(22) 


1.4(22) 


1.4(22) 


1.4(22) 


1.0(22) 


H + 


5.7(15) 


7.9(15) 


2.8(16) 


2.0(17) 


1.3(18) 


H 2 + 


2.7(09) 


2.1(10) 


2.0(11) 


2.0(12) 


1.6(13) 


H+ 


4.0(12) 


4.2(13) 


4.0(14) 


3.2(15) 


6.7(15) 


e 


8.2(17) 


8.2(17) 


8.5(17) 


1.0(18) 


2.3(18) 


CH + 


3.2(12) 


3.2(12) 


3.1(12) 


2.9(12) 


2.3(12) 


CH+ 


1.0(12) 


1.0(12) 


1.0(12) 


1.0(12) 


9.8(11) 


CH 


3.1(13) 


3.3(13) 


7.0(13) 


1.6(14) 


5.9(13) 


O 


6.6(18) 


6.6(18) 


6.7(18) 


7.0(18) 


8.1(18) 


+ 


1.7(12) 


2.4(12) 


8.6(12) 


5.8(13) 


3.5(14) 


OH + 


1.9(12) 


1.9(12) 


2.5(12) 


6.9(12) 


4.3(13) 


H 2 + 


6.2(11) 


6.5(11) 


9.3(11) 


3.6(12) 


4.2(13) 


H s O + 


1.7(14) 


2.5(14) 


3.5(14) 


6.6(14) 


6.1(14) 


OH 


7.2(14) 


1.7(15) 


8.6(15) 


5.1(16) 


1.7(17) 


H 2 


1.7(17) 


2.6(16) 


4.9(16) 


8.7(16) 


5.8(16) 


C + 


7.7(17) 


7.7(17) 


7.7(17) 


7.8(17) 


9.1(17) 


C 


4.5(17) 


4.8(17) 


5.3(17) 


6.7(17) 


1.0(18) 


CO 


3.2(18) 


3.2(18) 


3.2(18) 


3.0(18) 


2.5(18) 


CO + 


7.4(10) 


7.4(10) 


7.6(10) 


1.1(11) 


1.0(12) 


HCO + 


6.3(13) 


8.9(14) 


2.5(15) 


4.7(15) 


1.4(15) 


N + 


4.3(10) 


4.4(11) 


4.6(12) 


5.3(13) 


4.0(14) 


NH + 


3.1(11) 


3.1(11) 


2.9(11) 


2.5(11) 


3.8(11) 


NH+ 


1.1(11) 


1.1(11) 


1.2(11) 


2.0(11) 


7.3(11) 


HCN 


1.2(14) 


2.8(14) 


8.0(14) 


9.5(14) 


9.0(13) 


HNC 


7.5(13) 


1.1(14) 


2.5(14) 


3.1(14) 


8.9(13) 


CS 


8.6(11) 


5.2(12) 


2.3(13) 


5.0(13) 


1.0(13) 



Molecular hydrogen column densities only change significantly 
for the highest two CR rates, ~ 20 and ~ 40 percent lower, re- 
spectively. 

The H + column density is lower by a factor ~ 3 for the lowest 
and ~ 5 to ~ 15 percent higher for the highest two CR rates. The 
H + column density is dominated by + H — > H + + H2 in the 
radical region at f = 5 ■ 10~ 17 s _! and thus lower, while at higher 
CR rates the column density is increased because of the only 
partial transition to H2 

+ , OH, OH + , H 2 and H 2 + , and H 3 + : The trend seen 
for the + column density is very close to what is seen for H + , 
i.e, a factor of ~ 3 lower for ( - 5 ■ 10~ 17 s _1 , while slightly 
higher for the highest CR rates. Similar trends are seen for OH + 
andH 2 + . 

H2O is only significantly affected (lowered by a factor ~ 2) 
at the highest CR rate. OH column density is lower by a factor 
~ 3 for the lowest CR rates, and dominated by the abundance in 
the unshielded region of the cloud. This difference is much less 
(5-30 percent) for the higher CR rates. H 2 and H 3 + column 
densities are only affected for the higest two CR rates. 

In all we conclude that the overall trends are not affected by 
uncertainties in the H2 abundance, the main driver of the species 
of interest in our paper. The largest changes are seen when the 
integrated column density of a species is dominated by its abun- 
dance in the radical region, i.e., in the models with small CR 
rates. 



10 



R. Meijerink et al.: Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating. 



Table 5. Column densities for Model set 2 



Species Cosmic ray rate (s ') 





5(-17) 


5(-16) 


5(-15) 


5(-14) 


5(-13) 




Column densities (cm" 


2 ) 






H 


4.4(21) 


6.6(21) 


1.6(22) 


2.7 


22) 


3.2(22) 


H 2 


1.4(22) 


1.3(22) 


8.6(21) 


2.8 


21) 


5.2(20) 


H + 


5.6(16) 


4.2(17) 


3.3(18) 


2.8 


19) 


1.5(20) 


H 2 + 


6.3(11) 


6.0(12) 


4.5(13) 


1.8 


14) 


3.2(14) 


Ht 


1.1(15) 


5.8(15) 


1.9(15) 


3.1 


14) 


2.0(13) 


e 


9.5(17) 


1.4(18) 


5.5(18) 


3.4 


19) 


1.7(20) 


CH + 


7.9(10) 


4.1(11) 


5.8(11) 


1.3 


11) 


5.7(10) 


CH+ 


1.4(11) 


3.9(11) 


6.5(11) 


1.3 


11) 


6.1(09) 


CH 


4.5(13) 


1.7(14) 


6.0(13) 


4.3 


12) 


4.7(11) 


O 


6.8(18) 


7.8(18) 


1.1(19) 


1.1 


19) 


1.1(19) 


+ 


1.3(13) 


1.1(14) 


6.8(14) 


3.8 


15) 


1.7(16) 


OH + 


1.2(12) 


9.8(12) 


5.7(13) 


3.0 


14) 


2.5(14) 


H 2 + 


7.3(11) 


6.5(12) 


2.3(13) 


3.8 


13) 


1.6(12) 


H,0 + 


5.7(14) 


4.3(14) 


1.4(13) 


1.8 


12) 


3.0(09) 


OH 


1.8(16) 


5.9(16) 


7.5(15) 


3.3 


15) 


1.0(14) 


H 2 


1.1(17) 


1.5(17) 


1.6(15) 


1.7 


14) 


2.7(11) 


C + 


8.4(17) 


8.7(17) 


1.8(18) 


3.5 


18) 


4.1(18) 


C 


6.8(17) 


9.3(17) 


2.6(18) 


9.9 


17) 


3.6(17) 


CO 


2.9(18) 


2.7(18) 


4.0(16) 


1.8 


15) 


3.9(12) 


co + 


7.8(09) 


7.3(10) 


2.9(11) 


2.4 


11) 


4.3(09) 


HCO + 


8.2(14) 


6.0(14) 


5.0(12) 


1.7 


11) 


9.7(07) 


N + 


9.5(12) 


4.8(13) 


3.4(14) 


2.0 


15) 


6.9(15) 


NH + 


1.8(11) 


1.8(11) 


2.2(11) 


1.3 


12) 


4.8(12) 


NH+ 


2.6(10) 


1.6(11) 


2.3(11) 


3.0 


11) 


4.9(10) 


HCN 


2.2(14) 


1.2(14) 


2.6(12) 


1.5 


10) 


3.5(06) 


HNC 


2.5(14) 


1.8(14) 


6.1(12) 


3.4 


10) 


3.3(06) 


CS 


1.8(14) 


6.2(13) 


2.5(10) 


4.9 


06) 


3.2(03) 



4. Integrated line fluxes 

Line fluxes for a number of species discussed in the previous 
section are shown in Tables [6] through [9] Here we highlight the 
main results: 

Fine-structure lines: The two important cooling lines, 
[OI] 63 /mi and [CII] 158 /mi, that mainly determine the energy 
budget of interstellar clouds, have very similar fluxes for differ- 
ent CR rates in the high density case. [OI] 63 //m is the same 
within a factor of 2 in the high density case for all mechanical 
heating and CR rates, but does increase a factor 5 in the low den- 
sity case. The critical density for excitation of the [OI] 63 /mi is 
Kent — 5 • 10 5 ctrT 3 , and the low density model is not even close 
to being thermalized, and temperature differences cause a large 
change in emitted flux. The [CII] 158 //m has a very low critical 
density n cr i, = 2.8 ■ 10 3 cm -3 , and the increase in fluxes is not 
explained by the rise in temperature at small column densities. 
Additional flux is emitted in low density clouds, because there 
is no full transition to neutral carbon and CO, for the highest 
CR rates (see Section 2). This also explains why the two neutral 
carbon lines, [CI] 610 //m and [CI] 370 //m, show an increase in 
the high density models. The low density models have a max- 
imum intensity at £ ~ 5 ■ 10~ 15 s _1 , since for higher CR rates 
the main carbon budget remains in C + . The two ionized nitrogen 
lines, [Nil] 205 /mi and [Nil] 122 /mi, increase by 4 (3) orders of 
magnitude in the high (low) density case. Especially low density 
clouds are expected to yield detectable fluxes: the [Nil] 122 /mi 
lines has only a 10 times weaker flux than the [CII] line at the 
highest CR rate. Unfortunately, + does not have fine-structure 
transitions in the mid and far-infrared. Available transitions are 
in the visible and UV, and have very high excitation energies. 

H^ 95 um: This transi tion is calcalated following the discus- 
sion in iPan & Okal {[986). The fluxes in the high density case 




Column density [cm 2 ] 



Fig. 8. Abundances of hydrogen and water chemistry for an 
XDR model with density n - 10 5 5 crrr 3 and F x - 
28.5 erg cirT 2 s _1 . 



are very similar to the [Nil] 122 /mi, not extremely strong, but 
potentially detectable, when clouds are exposed to very large CR 
fluxes. 

HCN, HNC, and HCO + : The HCN and HNC do not seem 
very helpful in the discussion of cosmic ray rates, as there are 
no obvious trends. The HCN/HNC ratio does sh ow a response 
to me chanical heating as already pointed out in lLoenen et al.l 
(2008), and especially HCN is boosted by a few orders of mag- 
nitude. An exception is the high density with the high CR rate, 
where He + destroys the HCN and HNC very effectively. Also 
HCO + does not exhibit very obvious trends, except for very high 
mechanical heating rates (model set lc). 

CO and l3 CO: In the high density models, there is very lit- 
tle response when increasing the cosmic ray flux. Fluxes are 
boosted, when mechanical heating plays an important role in the 
energy budget of the gas, but, contrary to the HCN and HNC 
lines, it only affects high-/ transitions. When 13 CO line fluxes 
are also available, mechanical heating effects are potentially also 
seen in the lower-/ transitions (/ < 3). The CO to 13 CO line ra- 
tio increases, indicating smaller optical depths, while CO fluxes 
only decrease by a factor 2 



5. Key Diagnostics 

The aim of this paper is to pinpoint potential line diagnostics that 
can be used to trace an enhanced cosmic ray radiation field and 
mechanical heating in the presence of starformation such as the 
central regions of Arp 220 and NGC 253. One of the remark- 
able results is that many commonly observed lines, such as, the 
[OI] 63 /mi, [CI] 609 /mi, the [CII] 158 /mi and the low-/ CO 
lines, show very moderate responses to high cosmic ray rates 
and mechanical feedback, when the medium is already exposed 
to a large amount of UV. This is a very useful result, since these 
lines can serve as a reference point to other lines fluxes, that are 
impacted. 



R. Meijerink et al.: Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating. 
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Line 



Cosmic ray rate (s ') 



Line 



Cosmic ray rate (s ') 







5(-17) 


5(- 


16) 


5(-15) 


5(-14) 


5(-13) 




5(-17) 


5(-16) 


5(-15) 


5(-14) 


5(-13) 






Line intensities (er 


gs" J 


cm 2 


sr 1 ) 












Line intensities (erg s 


cm 2 


sr 1 ) 










[CII] 158 fim 




l.l(-3) 


l.K 


-3) 


l.K 


-3) 


l.K 


-3) 


1.3( 


-3) 


[CII] 158 yum 


l.K-3) 


l.l(-3) 


l.K 


-3) 


1.2( 


-3) 


1.3( 


-3) 


[CI] 610 yum 




3.7(-6) 


3.9( 


-6) 


4.2( 


-6) 


4.8( 


-6) 


l.K 


-5) 


[CI] 610 fim 


3.4(-6) 


3.6(-6) 


3.9( 


-6) 


4.9( 


-6) 


7.3( 


-6) 


[CI] 370 yum 




2.1(-5) 


2.2( 


-5) 


2.3( 


-5) 


2.6( 


-5) 


6.2( 


-5) 


[CI] 370 urn 


2.4(-5) 


2.5(-5) 


2.8( 


-5) 


3.4( 


-5) 


5.2( 


-5) 


[OI] 63 fim 




7.8(-2) 


7.8( 


-2) 


7.8( 


-2) 


8.2( 


-2) 


9.2( 


-2) 


[OI] 63 fim 


9.3(-2) 


9.4(-2) 


9.4( 


-2) 


9.8( 


-2) 


l.K 


-1) 


[OI] 146 fim 




2.3(-3) 


2.3( 


-3) 


2.3( 


-3) 


2.4( 


-3) 


2.9( 


-3) 


[OI] 146 yum 


3.7(-3) 


3.7(-3) 


3.8( 


-3) 


3.9( 


-3) 


4.1( 


-3) 


[Nil] 205 iim 




4.4(-ll) 


4.4( 


10) 


4.7( 


-9) 


5.5( 


-8) 


4.2( 


-7) 


[Nil] 205 yum 


4.7(-ll) 


4.7(-10) 


5.0( 


-9) 


5.7( 


-8) 


4.3( 


-7) 


[Nil] 122 yum 




6.6(-10) 


6.6( 


-9) 


7.K 


-8) 


8.4( 


-7) 


6.6( 


-6) 


[Nil] 122 yum 


7.2(-10) 


7.2(-9) 


7.7( 


-8) 


8.8( 


-7) 


6.7( 


-6) 


H+ 95 yum 




7.0(-10) 


6.2( 


-9) 


6.0( 


-8) 


6.51 


-7) 


1.9( 


-6) 


H+ 95 yum 


3.6(-9) 


3.8(-8) 


3.6( 


-7) 


2.8( 


-6) 


6.0( 


-6) 


HCN 7 = 1 - 





1.5(-8) 


1.8( 


-8) 


1.8( 


-8) 


2.9( 


-8) 


l.K 


-8) 


HCN 7=1-0 


3.1(-8) 


5.8(-8) 


1.2( 


-7) 


1.3( 


-7) 


2.0( 


-8) 


HCN 7 = 2- 


1 


1.3(-7) 


1.5( 


-7) 


1.5( 


-7) 


2.2( 


-7) 


1.3( 


-7) 


HCN 7 = 2-1 


2.9(-7) 


4.3(-7) 


6.7( 


-7) 


7.4( 


-7) 


2.8( 


-7) 


HCN 7 = 3- 


2 


3.2(-7) 


3.7( 


-7) 


3.6( 


-7) 


5.81 


-7) 


3.3( 


-7) 


HCN 7 = 3-2 


8.3(-7) 


1.2(-6) 


1.9( 


-6) 


2.0( 


-6) 


8.K 


-7) 


HCN 7 = 4- 


3 


3.8(-7) 


4.7( 


-7) 


4.5( 


-7) 


8.9( 


-7) 


4.4( 


-7) 


HCN 7 = 4-3 


1.6(-6) 


2.6(-6) 


4.1( 


-6) 


4.4( 


-6) 


1.5( 


-6) 


HCN 7 = 5- 


4 


3.3(-7) 


4.0( 


-7) 


3.9( 


-7) 


7.7( 


-7) 


4.0( 


-7) 


HCN 7 = 5-4 


2.0(-6) 


4.2(-6) 


7.4( 


-6) 


8.K 


-6) 


2.0( 


-6) 


HCN 7 = 6- 


5 


2.6(-7) 


3.0( 


-7) 


3.0( 


-7) 


5.2( 


-7) 


3.K 


-7) 


HCN 7 = 6-5 


1.8(-6) 


5.2(-6) 


1.2( 


-5) 


1.3( 


-5) 


1.9( 


-6) 


HCN 7 = 7- 


6 


2.0(-7) 


2.2( 


-7) 


2.2( 


-7) 


3.2( 


-7) 


2.2( 


-7) 


HCN 7 = 7-6 


1.5(-6) 


4.6(-6) 


1.5( 


-5) 


1.8( 


-5) 


1.6( 


-6) 


HNC 7=1- 





Ll(-8) 


1.3( 


-8) 


2.0( 


-8) 


3.6( 


-8) 


1.7( 


-8) 


HNC 7=1-0 


2.2(-8) 


3.0(-8) 


5.5( 


-8) 


6.0( 


-8) 


2.2( 


-8) 


HNC 7 = 2- 


1 


9.9(-8) 


l.K 


-7) 


1.6( 


-7) 


2.5( 


-7) 


1.7( 


-7) 


HNC 7 = 2-1 


2.3(-7) 


2.9(-7) 


4.3( 


-7) 


4.7( 


-7) 


2.8( 


-7) 


HNC 7 = 3- 


2 


2.1(-7) 


2.4( 


-7) 


3.9( 


-7) 


6.7( 


-7) 


4.5( 


-7) 


HNC 7 = 3-2 


6.5(-7) 


8.2(-7) 


1.2( 


-6) 


1.3( 


-6) 


8.K 


-7) 


HNC 7 = 4- 


3 


2.3(-7) 


2.6( 


-7) 


4.9( 


-7) 


l.K 


-6) 


6.3( 


-7) 


HNC 7 = 4-3 


l.K-6) 


1.5(-6) 


2.5( 


-6) 


2.8( 


-6) 


1.5( 


-6) 


HNC 7 = 5- 


4 


1.9(-7) 


2.0( 


-7) 


3.4( 


-7) 


l.K 


-6) 


5.7( 


-7) 


HNC 7 = 5-4 


L2(-6) 


2.0(-6) 


4.0( 


-6) 


4.5( 


-6) 


1.9( 


-6) 


HNC 7 = 6- 


5 


1.4(-7) 


1.5( 


-7) 


2.2( 


-7) 


6.K 


-7) 


4.2( 


-7) 


HNC 7 = 6-5 


l.K-6) 


1.8(-6) 


4.7( 


-6) 


5.6( 


-6) 


1.8( 


-6) 


HNC 7 = 7- 


6 


9.7(-8) 


1.0( 


-7) 


1.3( 


-7) 


3.3( 


-7) 


2.8( 


-7) 


HNC 7 = 7-6 


9.2(-7) 


1.5(-6) 


3.8( 


-6) 


4.9( 


-6) 


1.5( 


-6) 


HCO + 7 = 1 


-0 


2.7(-8) 


1.0( 


-7) 


1.5( 


-7) 


1.6( 


-7) 


5.5( 


-8) 


HCO + 7=1-0 


l-9(-8) 


1.4(-7) 


3.3( 


-7) 


5.5( 


-7) 


1.9( 


-7) 


HCO + 7 = 2 


- 1 


2.8(-7) 


7.7( 


-7) 


9.8( 


-7) 


1.2( 


-6) 


7.3( 


-7) 


HCO+ 7 = 2-1 


3.8(-7) 


1.7(-6) 


2.9( 


-6) 


3.9( 


-6) 


2.5( 


-6) 


HCO + 7 = 3 


-2 


8.K-7) 


2.2( 


-6) 


2.7( 


-6) 


3.3( 


-6) 


3.K 


-6) 


HCO + 7 = 3-2 


1.3(-6) 


4.7(-6) 


6.5( 


-6) 


8.2( 


-6) 


6.6( 


-6) 


HCO+ 7 = 4 


-3 


L4(-6) 


4.4( 


-6) 


5.5( 


-6) 


6.8( 


-6) 


6.8( 


-6) 


HCO + 7 = 4-3 


2.7(-6) 


9.5(-6) 


1.3( 


-5) 


1.5( 


-5) 


1.3( 


-5) 


HCO + 7 = 5 


-4 


1.6(-6) 


7.2( 


-6) 


9.4( 


-6) 


1.2( 


-5) 


9.4( 


-6) 


HCO + 7 = 5-4 


4.K-6) 


1.7(-5) 


2.3( 


-5) 


2.7( 


-5) 


2.3( 


-5) 


HCO+ 7 = 6 


-5 


1.2(-6) 


1.0( 


-5) 


1.4( 


-5) 


1.9( 


-5) 


1.3( 


-5) 


HCO+ 7 = 6-5 


4.9(-6) 


2.7(-5) 


3.8( 


-5) 


4.5( 


-5) 


3.7( 


-5) 


HCO + 7 = 7 


-6 


8.2(-7) 


1.2( 


-5) 


2.0( 


-5) 


2.8( 


-5) 


1.5( 


-5) 


HCO + 7 = 7-6 


4.5(-6) 


3.9(-5) 


5.9( 


-5) 


7.K 


-5) 


5.5( 


-5) 


CO 7 = 1 - 




3.0(-7) 


3.0( 


-7) 


3.K 


-7) 


3.2( 


-7) 


3.K 


-7) 


CO 7 = 1 - 


4.2(-7) 


4.1(-7) 


4.1( 


-7) 


4.0( 


-7) 


3.6( 


-7) 


CO 7 = 2 - 1 




2.9(-6) 


3.0( 


-6) 


3.0( 


-6) 


3.3( 


-6) 


3.9( 


-6) 


CO 7 = 2 - 1 


6.1(-6) 


6.1(-6) 


6.2( 


-6) 


6.K 


-6) 


5.8( 


-6) 


CO 7 = 3 - 2 




1.0(-5) 


1.0( 


-5) 


l.K 


-5) 


1.2( 


-5) 


1.5( 


-5) 


CO 7 = 3 - 2 


2.5(-5) 


2.5(-5) 


2.51 


-5) 


2.5( 


-5) 


2.4( 


-5) 


CO 7 = 4 - 3 




2.4(-5) 


2.4( 


-5) 


2.5( 


-5) 


2.8( 


-5) 


3.9( 


-5) 


CO 7 = 4 - 3 


6.3(-5) 


6.3(-5) 


6.4( 


-5) 


6.4( 


-5) 


6.4( 


-5) 


CO 7 = 5 - 4 




4.3(-5) 


4.4( 


-5) 


4.6( 


-5) 


5.3( 


-5) 


7.8( 


-5) 


CO 7 = 5 - 4 


1.3(-4) 


1.3(-4) 


1.3( 


-4) 


1.3( 


-4) 


1.3( 


-4) 


CO 7 = 6 - 5 




6.7(-5) 


6.8( 


-5) 


7.K 


-5) 


8.4( 


-5) 


1.3( 


-4) 


CO 7 = 6 - 5 


2.2(-4) 


2.2(-4) 


2.3( 


-4) 


2.3( 


-4) 


2.3( 


-4) 


CO 7 = 7 - 6 




9.2(-5) 


9.4( 


-5) 


9.8( 


-5) 


1.2( 


-4) 


2.0( 


-4) 


CO 7 = 7 - 6 


3.5(-4) 


3.6(-4) 


3.6( 


-4) 


3.6( 


-4) 


3.71 


-4) 


CO 7 = 8 - 7 




1.1 (-4) 


l.K 


-4) 


1.2( 


-4) 


1.5( 


-4) 


2.8( 


-4) 


CO 7 = 8 - 7 


5.2(-4) 


5.2(-4) 


5.3( 


-4) 


5.3( 


-4) 


5.5( 


-4) 


CO 7 = 9 - 8 




1.2(-4) 


1.3( 


-4) 


1.3( 


-4) 


1.7( 


-4) 


3.6( 


-4) 


CO 7 = 9 - 8 


7.K-4) 


7.3(-4) 


7.3( 


-4) 


7.3( 


-4) 


7.7( 


-4) 


CO 7 = 10 - 


9 


1.1 (-4) 


1.2( 


-4) 


1.2( 


-4) 


1.8( 


-4) 


4.3( 


-4) 


CO 7 = 10-9 


9.4(-4) 


9.6(-4) 


9.7( 


-4) 


9.7( 


-4) 


1.0( 


-3) 


CO 7= 11 - 


10 


8.2(-5) 


8.9( 


-5) 


9.4( 


-5) 


1.6( 


-4) 


4.8( 


-4) 


C07= 11 - 10 


1.2(-3) 


1.2(-3) 


1.2( 


-3) 


1.2( 


-3) 


1.3( 


-3) 


CO 7 = 12 - 


11 


3.8(-5) 


4.3( 


-5) 


4.8( 


-5) 


l.K 


-4) 


5.0( 


-4) 


CO 7= 12- 11 


1.5(-3) 


1.5(-3) 


1.5( 


-3) 


1.5( 


-3) 


1.6( 


-3) 


CO 7 = 13 - 


12 


1.2(-5) 


1.4( 


-5) 


1.6( 


-5) 


5.K 


-5) 


4.7( 


-4) 


CO 7 = 13 - 12 


1.7(-3) 


1.8(-3) 


1.81 


-3) 


1.7( 


-3) 


1.8( 


-3) 


CO 7 = 14 - 


13 


4.7(-6) 


5.2( 


-6) 


5.7( 


-6) 


1.8( 


-5) 


3.8( 


-4) 


CO 7 = 14- 13 


2.0(-3) 


2.0(-3) 


2.0< 


-3) 


2.0( 


-3) 


2.K 


-3) 


CO 7 = 17 - 


16 


3.4(-6) 


3.4( 


-6) 


3.4( 


-6) 


3.8( 


-6) 


5.8( 


-5) 


CO 7 = 17 - 16 


2.2(-3) 


2.4(-3) 


2.3( 


-3) 


2.2( 


-3) 


2.3( 


-3) 


CO 7 = 20 - 


19 


3.6(-6) 


3.6( 


-6) 


3.6( 


-6) 


3.5( 


-6) 


7.8( 


-6) 


CO 7 = 20 - 19 


9.K-4) 


l.l(-3) 


1.0( 


-3) 


8.7( 


-4) 


1.2( 


-3) 


CO 7 = 24 - 


23 


2.4(-6) 


2.4( 


-6) 


2.3( 


-6) 


2.2( 


-6) 


2.4( 


-6) 


CO 7 = 24 - 23 


2.7(-5) 


3.7(-5) 


3.K 


-5) 


2.6( 


-5) 


5.51 


-5) 


CO 7 = 30 - 


29 


4.9(-7) 


4.9( 


-7) 


4.8( 


-7) 


4.4( 


-7) 


4.9( 


-7) 


CO 7 = 30 - 29 


8.0(-7) 


8.1(-7) 


7.81 


-7) 


7.2( 


-7) 


7.8( 


-7) 


CO 7 = 35 - 


34 


1.0(-7) 


1.0( 


-7) 


1.0( 


-7) 


9.9( 


-8) 


1.3( 


-7) 


CU J = 33 — 34 


1.7(-7) 


1.7(-7) 


1.7( 


-7) 


1.6( 


-7) 


1.9( 


-7) 


13 CO 7=1- 





4.1(-8) 


4.1( 


-8) 


4.0( 


-8) 


3.6( 


-8) 


1.4( 


-8) 


CU 7 = 1—0 


1.5(-8) 


1.4(-8) 


1.4( 


-8) 


1.4( 


-8) 


l.K 


-8) 


13 CO 7 = 2- 


1 


7.9(-7) 


7.9( 


-7) 


7.9( 


-7) 


7.6( 


-7) 


3.7( 


-7) 


CO 7 = 2-1 


4.3(-7) 


4.2(-7) 


4.2( 


-7) 


4.0( 


-7) 


3.2( 


-7) 


13 CO 7 = 3- 


2 


3.4(-6) 


3.4( 


-6) 


3.4( 


-6) 


3.5( 


-6) 


2.2( 


-6) 


CO 7 = 3-2 


2.9(-6) 


2.8(-6) 


2.8( 


-6) 


2.7( 


-6) 


2.2( 


-6) 


13 CO 7 = 4- 


3 


8.4(-6) 


8.5( 


-6) 


8.5( 


-6) 


9.1( 


-6) 


7.0( 


-6) 


CO 7 = 4-3 


1.0(-5) 


1.0(-5) 


1.0( 


-5) 


9.6( 


-6) 


7.9( 


-6) 


CO 7 = 5- 


4 


1.5(-5) 


1.5( 


-5) 


1.5( 


-5) 


1.7 


-5) 


1.5( 


-5) 


CO 7 = 5-4 


2.6(-5) 


2.5(-5) 


2.5( 


-5) 


2.4( 


-5) 


2.0( 


-5) 


\ A J J — U 


j 


2.2(-5) 


2.2( 


j ) 


2.2( 


-5) 


2.6( 


-5) 


2.7( 


-5) 




5.1(-5) 


5.0(-5) 


5.0( 


-5) 


4.8( 


-5) 


4.K 


-5) 


13 CO 7 = 7- 


6 


2.5(-5) 


2.6( 


-5) 


2.6( 


-5) 


3.3( 


-5) 


3.9( 


-5) 


13 CO 7 = 7-6 


8.5(-5) 


8.5(-5) 


8.5( 


-5) 


8.K 


-5) 


7.0( 


-5) 


13 C07 = 8- 


7 


2.3(-5) 


2.4( 


-5) 


2.5( 


-5) 


3.4( 


-5) 


5.0( 


-5) 


13 C07 = 8-7 


1.3(-4) 


1.3(-4) 


1.3( 


-4) 


1.2( 


-4) 


l.K 


-4) 


13 C07 = 9- 


8 


1.6(-5) 


1.7( 


-5) 


1.8( 


-5) 


2.8( 


-5) 


5.6( 


-5) 


13 C0 7 = 9-8 


L7(-4) 


1.7(-4) 


1.71 


-4) 


1.7( 


-4) 


1.5( 


-4) 


13 CO 7 = 10 


-9 


8.3(-6) 


9.K 


-6) 


9.7( 


-6) 


1.8( 


-5) 


5.5( 


-5) 


13 C0 7 = 10-9 


2.K-4) 


2.1 (-4) 


2.K 


-4) 


2.0( 


-4) 


1.9( 


-4) 


13 C0 7= 11 


- 10 


3.4(-6) 


3.8( 


-6) 


4.2( 


-6) 


9.6( 


-6) 


4.7( 


-5) 


13 C0 7 = 11 - 10 


2.4(-4) 


2.5(-4) 


2.4( 


-4) 


2.3( 


-4) 


2.2( 


-4) 


13 CO 7 = 12 


- 11 


L2(-6) 


1.4( 


-6) 


1.5( 


-6) 


4.2( 


-6) 


3.6( 


-5) 


13 C0 7 = 12-11 


2.504) 


2.6(-4) 


2.5( 


-4) 


2.4( 


-4) 


2.3( 


-4) 


13 C0 7 = 13 


- 12 


4.0(-7) 


4.4( 


-7) 


5.0( 


-7) 


1.6( 


-6) 


2.4( 


-5) 


13 C0 7 = 13 - 12 


2.4(-4) 


2.5(-4) 


2.4( 


-4) 


2.3( 


-4) 


2.4( 


-4) 


13 C0 7 = 14 


- 13 


1.5(-7) 


1.6( 


-7) 


1.8( 


-7) 


5.71 


-7) 


1.4< 


-5) 


13 C0 7= 14-13 


2.0(-4) 


2.2(-4) 


2.K 


-4) 


2.0( 


-4) 


2.2( 


-4) 
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Table 8. Line intensities for Model set lc Table 9. Line intensities for Model set 2 



Line 


Cosmic ray rate (s ') Line 
5(-17) 5(-16) 5(-15) 5(-14) 5(-13) 


Cosmic ray rate (s 
5(-17) 5(-16) 5(-15) 


- 1 ) 

5(-14) 


5(-13) 




Line intensities (erg s~ J crrr 2 sr -1 ) 


Line intensities (erg s~ J cirr 2 sr~') 







[CII] 158 fim 




1.1( 


-3) 


1.1( 


-3) 


1.1( 


-3) 


1.1( 


-3) 


1.3( 


-3) 


[CI] 610 yum 




8.3( 


-7) 


1.3( 


-6) 


2.5( 


-6) 


3.8( 


-6) 


7.7( 


-6) 


[CI] 370 fim 




7.1( 


-6) 


1.1( 


-5) 


2.1( 


-5) 


3.2( 


-5) 


6.3( 


-5) 


[OI] 63 yum 




1.1( 


-1) 


1.1( 


-1) 


1.1( 


-1) 


1.3( 


-1) 


1.5( 


-1) 


[OI] 146 fim 




3.9( 


-3) 


3.9( 


-3) 


4.1( 


-3) 


5.4( 


-3) 


6.7( 


-3) 


[Nil] 205 yum 




7.3(- 


11) 


7.3(- 


10) 


7.5( 


-9) 


7.1( 


-8) 


4.7( 


-7) 


[Nil] 122 fim 




1.1( 


-9) 


1.1( 


-8) 


1.1( 


-7) 


1.1( 


-6) 


7.4( 


-6) 


H+ 95 yum 




1.8( 


-9) 


1.7( 


-8) 


1.8( 


-7) 


2.1( 


-6) 


9.9( 


-6) 


HCN 7 = 1 - 





3.9( 


-5) 


4.0( 


-5) 


2.2( 


-5) 


4.7( 


-6) 


1.8( 


-7) 


HCN 7 = 2- 


1 


2.1( 


-5) 


2.0( 


-5) 


1.7( 


-5) 


9.5( 


-6) 


2.7( 


-6) 


HCN 7 = 3 - 


2 


1.8( 


-5) 


1.9( 


-5) 


1.8( 


-5) 


1.4( 


-5) 


4.1( 


-6) 


HCN 7 = 4 - 


3 


3.0( 


-5) 


3.1( 


-5) 


3.0( 


-5) 


2.6( 


-5) 


1.2( 


-5) 


HCN 7 = 5 - 


4 


5.2( 


-5) 


5.3( 


-5) 


5.3( 


-5) 


4.4( 


-5) 


2.0( 


-5) 


HCN 7 = 6- 


5 


9.1( 


-5) 


9.3( 


-5) 


9.3( 


-5) 


7.5( 


-5) 


3.5( 


-5) 


HCN 7 = 7 - 


6 


1.5( 


-4) 


1.6( 


-4) 


1.6( 


-4) 


1.2( 


-4) 


5.4( 


-5) 


HNC 7=1- 





7.9( 


-7) 


6.0( 


-7) 


4.4( 


-7) 


1.2( 


-7) 


L7( 


-8) 


HNC 7 = 2- 


1 


3.3( 


-6) 


2.9( 


-6) 


2.4( 


-6) 


1.1( 


-6) 


3.0( 


-7) 


HNC 7 = 3- 


2 


6.4( 


-6) 


5.9( 


-6) 


5.2( 


-6) 


4.4( 


-6) 


1.0( 


-6) 


HNC 7 = 4- 


3 


1.3( 


-5) 


1.2( 


-5) 


1.1( 


-5) 


7.4( 


-6) 


2.1( 


-6) 


HNC 7 = 5- 


4 


2.4( 


-5) 


2.2( 


-5) 


2.0( 


-5) 


1.3( 


-5) 


3.1( 


-6) 


HNC 7 = 6- 


5 


4.1( 


-5) 


3.8( 


-5) 


3.4( 


-5) 


2.0( 


-5) 


3.4( 


-6) 


HNC 7 = 7- 


6 


6.8( 


-5) 


6.2( 


-5) 


5.5( 


-5) 


3.0( 


-5) 


3.3( 


-6) 


HCO + 7 = 1 


-0 


4.9( 


-9) 


1.0( 


-8) 


3.9( 


-8) 


1.9( 


-7) 


4.9( 


-7) 


HCO + 7 = 2 


- 1 


1.4( 


-7) 


2.7( 


-7) 


8.8( 


-7) 


3.1( 


-6) 


7.6( 


-6) 


HCO + 7 = 3 


-2 


6.6( 


-7) 


1.2( 


-6) 


3.2( 


-6) 


8.7( 


-6) 


1.7( 


-5) 


HCO + 7 = 4 


-3 


1.7( 


-6) 


2.7( 


-6) 


6.8( 


-6) 


1.7( 


-5) 


3.1( 


-5) 


HCO + 7 = 5 


-4 


2.7( 


-6) 


4.2( 


-6) 


1.2( 


-5) 


3.0( 


-5) 


5.0( 


-5) 


HCO+ 7 = 6 


-5 


3.2( 


-6) 


5.1( 


-6) 


1.7( 


-5) 


4.7( 


-5) 


7.9( 


-5) 


HCO+ 7 = 7 


-6 


3.3( 


-6) 


5.2( 


-6) 


2.2( 


-5) 


7.0( 


-5) 


1.2( 


-4) 


CO 7 = 1 - 




2.1( 


-7) 


2.3( 


-7) 


2.7( 


-7) 


2.9( 


-7) 


2.3( 


-7) 


CO 7 = 2 - 1 




5.5( 


-6) 


5.9( 


-6) 


6.7( 


-6) 


6.9( 


-6) 


5.7( 


-6) 


CO 7 = 3 - 2 




3.1( 


-5) 


3.3( 


-5) 


3.6( 


-5) 


3.7( 


-5) 


3.2( 


-5) 


CO 7 = 4 - 3 




9.6( 


-5) 


1.0( 


-4) 


1.1( 


-4) 


1.1( 


-4) 


9.8( 


-5) 


CO 7 = 5 -4 




2.2( 


-4) 


2.2( 


-4) 


2.4( 


-4) 


2.4( 


-4) 


2.2( 


-4) 


CO 7 = 6 - 5 




4.1( 


-4) 


4.2( 


-4) 


4.4( 


-4) 


4.5( 


-4) 


4.2( 


-4) 


CO 7 = 7 - 6 




6.8( 


-4) 


7.0( 


-4) 


7.3( 


-4) 


7.5( 


-4) 


7.0( 


-4) 


CO 7 = 8-7 




1.0( 


-3) 


1.1( 


-3) 


1.1( 


-3) 


1.2( 


-3) 


1.1( 


-3) 


CO 7 = 9 - 8 




1.5( 


-3) 


1.5( 


-3) 


1.6( 


-3) 


1.6( 


-3) 


1.6( 


-3) 


CO 7 = 10- 


9 


2.0( 


-3) 


2.1( 


-3) 


2.2( 


-3) 


2.2( 


-3) 


2.1( 


-3) 


CO 7= 11 - 


10 


2.6( 


-3) 


2.7( 


-3) 


2.9( 


-3) 


2.9( 


-3) 


2.8( 


-3) 


CO 7= 12- 


11 


3.3( 


-3) 


3.4( 


-3) 


3.6( 


-3) 


3.7( 


-3) 


3.6( 


-3) 


C0 7= 13- 


12 


4.0( 


-3) 


4.1( 


-3) 


4.4( 


-3) 


4.5( 


-3) 


4.5( 


-3) 


CO 7 = 14 - 


13 


4.7( 


-3) 


4.9( 


-3) 


5.2( 


-3) 


5.4( 


-3) 


5.4( 


-3) 


CO 7= 17- 


16 


6.5( 


-3) 


6.8( 


-3) 


7.3( 


-3) 


7.6( 


-3) 


8.1( 


-3) 


CO 7 = 20 - 


19 


6.6( 


-3) 


7.0( 


-3) 


7.7( 


-3) 


8.3( 


-3) 


9.6( 


-3) 


CO 7 = 24 - 


23 


1.8( 


-3) 


2.1( 


-3) 


2.9( 


-3) 


3.7( 


-3) 


6.3( 


-3) 


CO 7 = 30 - 


29 


3.7( 


-5) 


4.2( 


-5) 


5.4( 


-5) 


6.9( 


-5) 


2.1( 


-4) 


CO 7 = 35 - 


34 


3.1( 


-6) 


3.2( 


-6) 


3.5( 


-6) 


3.6( 


-6) 


7.6( 


-6) 


13 CO 7=1- 





5.5( 


-9) 


6.1( 


-9) 


7.2( 


-9) 


7.6( 


-9) 


5.6( 


-9) 


13 CO 7 = 2- 


1 


1.7( 


-7) 


1.9( 


-8) 


2.3( 


-7) 


2.4( 


-7) 


1.7( 


-7) 


13 CO 7 = 3- 


2 


1.2( 


-6) 


1.4( 


-6) 


1.6( 


-6) 


1.7( 


-6) 


1.31 


-6) 


13 CO 7 = 4- 


3 


5.0( 


-6) 


5.5( 


-6) 


6.5( 


-6) 


6.8( 


-6) 


5.0( 


-6) 


13 CO 7 = 5- 


4 


1.4( 


-5) 


1.5( 


-5) 


1.8( 


-5) 


1.9( 


-5) 


1.4( 


-5) 


13 CO 7 = 6- 


5 


3.2( 


-5) 


3.5( 


-5) 


4.1( 


-5) 


4.3( 


-5) 


3.2( 


-5) 


13 CO 7 = 7- 


6 


6.1( 


-5) 


6.6( 


-5) 


7.8( 


-5) 


8.2( 


-5) 


6.3( 


-5) 


13 CO 7 = 8- 


7 


1.0( 


-4) 


1.1( 


-4) 


1.3( 


-4) 


1.4( 


-4) 


1.1( 


-4) 


13 CO 7 = 9- 


8 


1.6( 


-4) 


1.7( 


-4) 


2.0( 


-4) 


2.1( 


-4) 


1.7( 


-4) 


13 CO 7 = 10 


-9 


2.2( 


-4) 


2.4( 


-4) 


2.8( 


-4) 


3.0( 


-4) 


2.4( 


-4) 


13 C0 7 = 11 


- 10 


2.9( 


-4) 


3.1( 


-4) 


3.6( 


-4) 


3.9( 


-4) 


3.3( 


-4) 


13 CO 7 = 12 


- 11 


3.5( 


-4) 


3.8( 


-4) 


4.4( 


-4) 


4.8( 


-4) 


4.2( 


-4) 


13 CO 7 = 13 


- 12 


4.0( 


-4) 


4.3( 


-4) 


5.1( 


-4) 


5.5( 


-4) 


5.1( 


-4) 


13 C0 7 = 14 


- 13 


4.2( 


-4) 


4.6( 


-4) 


5.4( 


-4) 


6.0( 


-4) 


5.9( 


-4) 



[CII] 158 fim 




4.4(-4) 


4.7(-4) 


5.8(-4) 


8.6(-4) 


1.2(-3) 


[CI] 610 fim 




3.6(-6) 


4.6(-6) 


9.7(-6) 


6.6(-6) 


2.8(-6) 


[CI] 370 fim 




1.0(-5) 


1.3(-5) 


3.5(-5) 


3.0(-5) 


1.5(-5) 


[OI] 63 fim 




4.5(-4) 


6.4(-4) 


1.3(-3) 


2.4(-3) 


2.9(-3) 


[OI] 146 yum 




5.5(-6) 


7.0(-6) 


l.l(-5) 


2.K-5) 


3.1(-5) 


[Nil] 205 fim 




6.6(-9) 


5.3(-8) 


3.7(-7) 


2.1(-6) 


7.6(-6) 


[Nil] 122 fim 




4.2(-7) 


7.2(-7) 


5.5(-6) 


3.3(-5) 


1.1 (-4) 


H+ 95 yum 




4.7(-10) 


7.0(-9) 


3.2(-8) 


2.0(-8) 


3.6(-9) 


CO 7=1-0 




7.6(-8) 


1.3(-7) 


5.0(-8) 


2.8(-9) 


5.2(-12) 


CO 7 = 2 - 1 




4.0(-7) 


6.5(-7) 


4.4(-7) 


5.3(-8) 


1.2(-10) 


CO 7 = 3 - 2 




1.5(-6) 


2.0(-6) 


1.1 (-6) 


1.4(-7) 


3.6(-10) 


CO 7 = 4 - 3 




4.7(-6) 


5.7(-6) 


1.4(-6) 


1.4(-7) 


4.2(-10) 


CO 7 = 5 - 4 




1.2(-5) 


1.4(-5) 


7.6(-7) 


8.9(-8) 


2.9(-10) 


CO 7 = 6 - 5 




2.2(-5) 


2.6(-5) 


2.7(-7) 


4.1(-8) 


1.6(-10) 


CO 7 = 7 - 6 




3.2(-5) 


4.1(-5) 


7.5(-8) 


1.6(-8) 


7.0(-ll) 


CO 7 = 8 - 7 




3.4(-5) 


5.3(-5) 


1.9(-8) 


5.3(-9) 


2.8(-ll) 


CO 7 = 9 - 8 




2.3(-5) 


5.1(-5) 


4.4(-9) 


L6(-9) 


l.K-ll) 


CO 7 = 10 - 


9 


2.9(-6) 


2.7(-5) 


Ll(-9) 


4.9(-10) 


4.3(-12) 


CO 7 = 1 1 - 


10 


7.2(-9) 


1.8(-6) 


3.4(-10) 


1.5(-10) 


2.0(-12) 


CO 7 = 12 - 


1 1 


5.0(-12) 


4.0(-9) 


1.1(-10) 


5.2(-ll) 


1.1(-12) 


CO 7= 13- 


12 


3.6(-13) 


4.3(-ll) 


4.0(-ll) 


L9(-ll) 


6.5(-13) 


CO 7 = 14- 


13 


1.0(-13) 


2.1(-12) 


1.2(-11) 


5.8(-12) 


1.8(-13) 


13 CO 7=1- 





3.1(-8) 


4.6(-8) 


1.6(-9) 


5.K-11) 


9.1(-14) 


13 CO 7 = 2- 


1 


1.6(-7) 


2.6(-7) 


2.5(-8) 


1.2(-9) 


2.4(-12) 


13 CO 7 = 3- 


2 


3.3(-7) 


6.4(-7) 


5.3(-8) 


3.9(-9) 


9.5(-12) 


13 C0 7 = 4- 


3 


2.0(-7) 


7.9(-7) 


4.6(-8) 


5.0(-9) 


1.40-11) 


13 C07 = 5- 


4 


2.0(-8) 


2.5(-7) 


2.6(-8) 


3.9(-9) 


1-20-11) 


13 CO 7 = 6- 


5 


1.9(-9) 


4.4(-8) 


l.l(-8) 


2.3(-9) 


8.2(-12) 


13 CO 7 = 7- 


6 


1.2(-10) 


7.5(-9) 


4.4(-9) 


1.2(-9) 


4.7(-12) 


13 CO 7 = 8- 


7 


9.8(-12) 


LI (-9) 


1.5(-9) 


5.2(-10) 


2.4(-12) 



5.1. Tracers of Mechanical Heating 

Mechanical heating is only increasing the temperature of the 
gas, leaving the ionization fraction unaffected. Tracing this heat- 
ing mechanism with molecules that contain activation barriers 
in their formation route is thus very useful. Lines of, e.g., HCN 
and H2O are enhanced with respect to the aforementioned fine- 
structure lines, especially for the highest co nsidered mechanica l 
heating rate. It was already pointed out bv lLoenen et aD {2008) 
that the HCN/HNC line ratios increase when mechanical heating 
is important, since at high temperatures the HNC will be driven 
in to HCN (see also SectionEJ. Moreover, the HCN/CO line ra- 
tios also increase when mechanical heating effects are important, 
but the interpretation of an observed ratio is not straightforward, 
because the lines trace slightly different regions. 

5.2. Tracers of Enhanced Cosmic Ray Rates 

The [Nil] fine-structure lines show a very strong response to en- 
hanced cosmic ray rates, independent of mechanical heating ef- 
fects. The main concern with these lines, however, is that one has 
to disentangle the contributions arising from these CR exposed 
clouds and HII regions, that also contribute a significant amount 
of [Nil] emission in starforming environments. 

Most promising would be the study of lines from species im- 
portant in the water chemistry, especially when the ionized water 
related species are considered. The OH to H2O column density 
Tatio varies over two to three orders of magnitude for the cos- 
mic ray rates considered, from ~ 5 • 1CT 3 (^ = 5.0 • 10~ 17 s _1 ) 
to ~ 1 — 3 (£ = 5.0 • 10~ 13 s -1 ) in the high density case, and 
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~ 0.2 to ~ 400 in the low density case. So generally the OH lines 
are expected to be stronger in comparison to the water lines for 
higher cosmic ray rates. Fig. [8] shows an XDR model with den- 
sity n = 10 5 5 cirT 2 , and radiation field Fx = 28.5 erg ctrT 2 s" 1 
(a radiation field at ~170 pc distance from a typical AGN with 
Lx = 10 44 erg/s). From this we derive that the ratio varies be- 
tween ~ 8 (integrated up to A^h = 3 • 10 22 cm" 2 ) and ~ 1 
(A^h = 10 24 cm" 2 ). So for very high CR rates, the ratios become 
comparab le to a typic al X-ray illuminated cloud. 

IVan der Tak et alJ d2008l) already pointed out that the 
H3OVH2O ratio is affected by either cosmic rays or X-rays. 
The high density models show variations over 1 to 2 orders of 
magnitude in the high density case. Unfortunately, the interpre- 
tation is not entirely unambiguous, since the H^O + abundance 
decreases for the highest CR rates. Also, clouds illuminated by 
high CR rates have similar ratios as those exposed to X-rays, 
namely (1 - 6) • 10" 3 (£ = 5 • 10" 14 s" 1 ) and (5 - 10) • 10" 3 
(f = 5 ■ 10" 13 s" 1 ) for models set 1, and between 15 - 4 • 10" 3 
(depending on the size of the cloud) for the XDR model shown 
in Fig. [8] 

The current study reveals that additional observations of 
OH + and H20 + might help out in distinguishing very high CR 
rates from a typical XDR. The OH + and H20 + are significantly 
enhanced for CR rates that are 10 3 - 10 4 times larger than the 
accepted Milky Way value, but not as much as in the XDR. 
OH + /OH and H2OVH2O column density ratios can be as high 
as (3 - 5) ■ 10" 4 to (2 - 14) ■ 10" 4 , respectively, compared to 
(250-7.5) ■ 10" 4 and (64-4) ■ 10" 3 (N H = 3 ■ 10 22 - 10 24 cm" 2 ). 
XDRs and regions with enhanced CR rates can be distinguished 
when all species are considered. 

In this paper, no attempt is made t o model lines from these 
specie s. As already pointed out in, e.g.. iGonzalez- Alfonso et ail 
(2008), the radiation transfer of OH and H2O involves the con- 
sideration of both collisions and radiation pumping, which is be- 
yond the scope of this paper. Also, OH + and H20 + lines are not 
easily modeled. These species have formation and destruction 
times scales that are similar those for excitation. Besides that, 
collisional cross sections are unavailable at the moment, and 
quantum mechanical calculations are needed to obtain them. 
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6. Outlook 

OH + and H20 + have already been detected in a number of 
sources: 

- Ivan der Werf et all ((2010) show OH + and H20 + emission 
lines in the ULIRG Mrk231, which are about 30 percent of 
strength of the CO(5-6) to CO(13-12) lines. The paper show 
that these lines are the result from AGN activity in Mrk23 1 . 

- In the Milky Way, OH + has been detected in absorp- 
tion toward Sagittarius B2 using the A t acama Pathfinder 
Experiment (APEX) by IWvrowski et al.1 (1201 Ol). finding a 



column density of 2.4 ■ 10 15 cm" 2 . lOssenkopf et al.1 d2010l) 
found H 2 + in absorption toward, NGC 6334, DR21, again 
Sgr B2, obtaining column densities of 7.2 ■ 10 12 cm" 2 , 
2.3- 10 13 cm" 2 , and 1.1 - 10 15 cm" 2 , respectively. These galac- 
tic observations relate these species to the diffuse ISM. 

Extensive modeling of the excitation of water related species is 
very timely, since more and more observations will be available 
in a very short timefrace from observations with the recently 
lauched Herschel Space Telescope. 
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